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We address the quantum dot phase measurement problem in an open Aharonov-Bohm interfer- 
ometer, assuming multiple transport channels. In such a case, the quantum dot is characterized 
by more than one intrinsic phase for the electrons transmission. It is shown that the phase which 
would be extracted by the usual experimental method (i.e. by monitoring the shift of the Aharonov- 
Bohm oscillations, as in Schuster et al., Nature 385, 417 (1997)) does not coincide with any of the 
dot intrinsic phases, but is a combination of them. The formula of the measured phase is given. 
The particular case of a quantum dot containing a S = 1/2 spin is discussed and variations of the 
measured phase with less than -n are found, as a consequence of the multichannel transport. 



I. INTRODUCTION 



The recent developments in the building of mesoscopic devices allow to measure the phase (phase shift) of the 
electronic wave function after it is scattered by a quantum dot. Yacoby et al. [l[, placed a quantum dot in one of 
the arms of a quantum ring, realizing an Aharonov-Bohm (AB) interferometer and performed for the first time a 
phase measurement. In order to avoid multiple encirclements of the ring, Schuster et al. Q went a step further and 
opened the Aharonov-Bohm interferometer by adding a number of additional leads (see also 0, [E| and 0, 0] about 
the conditions of opening the interferometer). In this case, the electronic waves that reach the drain, have passed 
through the dot and through the reference arm only once, with the amplitudes tdot = \tdot\e l ^' d ° t and t re f = \t re f\e l ^ r " i 
respectively, yielding the total transmittancc through the interferometer: 

Tint = \tdot\ 2 + \tref\ 2 + 2 1 tdot | | t r ef \ COs((f>dot ~ <f>ref)- (1) 

The dot phase ipdot is extracted form Eq.l by varying the magnetic flux through the interferometer, and in 0] a 
phase lapse of 7r was found, between any pair of consecutive resonances. The expected behavior would have been a 
phase increase by n on each resonance (which Concordes with the experimental findings) , but also a constant value of 
the phase between resonances (when the dot occupancy does not change), as suggested by the Friedel sum rule 
The universal phase lapse behavior found by B is one of the longest-standing puzzles in mesoscopic physics, which 
lead to numerous theoretical studies -e.g. [41 lo-tlOl]- the issue being still under debate. 

In 2005, Avinun-Kalish et al. [llT ] launched a couple of new puzzles when performing phase measurements on a 
quantum dot that was initially empty and then gradually filled with electrons. Their experiment showed that the phase 
lapse of 7r was actually not universal, as found in but appears after the dot is occupied with several electrons. In 
the few-electrons regime, the phase variation shows a different behavior, also intriguing: one could notices a number 
of dips in the phase evolution of amplitude less than 7r, attributed by the authors to the Kondo effect. However, 
other experimental papers that addressed the phase measurement in the Kondo regime (e.g. (l2l |) reported a different 
behavior, namely an increase of the phase by 7r/2 on the resonance, which then remains constant at this value. The 



paper of Avinun-Kalish et al. [Tl[ opened a new trend of calculations on few-electrons systems (e.g. |13l - ll8| ). The 
reduced variation of the phase on and between (some) resonances, in the few-electrons regime, is a new puzzle in 
mesoscopic physics, and one of the motivations of our paper. We will show that the existence of multiple transport 
channels can lead to a reduced variation of the measured phase. 

To illustrate our statement, we shall address the scattering of an electron through a quantum dot with a S = 1/2 
magnetic impurity. The model has the advantage of being exactly solvable, while being a good approximation for the 
case when one electron is located in the dot, and another -itinerant- electron is scattered. Some authors proposed 
more complex models for impurity dots, e.g. [Tfj| - l23j . but this is not a purpose of this paper, where a simple model is 
needed to set the focus on the phase measurement procedure. 

The scattering of an electron on an impurity that possesses internal degrees of freedom was addressed before in 
mesoscopic physics, e.g by Mcllo et al. [24 



but not in the context of the interfcrometry and phase measurement. 
The outline of the paper is as follows: in section II wc describe the transport through the multichannel interferometer 
and derive the formula of the measured phase; sections III and IV address the case of a quantum dot with a magnetic 
impurity; the conclusions are given in section V. 
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II. THE MEASURED PHASE 



We consider an open Aharonov-Bohm interferometer, sketched in Fig.l. The term "open" is used in the sense of 
1 043 j meaning that multiple encirclements of the ring are neglected and only one interference is assumed between 
the partial electronic waves traveling through the upper and lower arms, respectively. The upper arm has an embedded 
quantum dot and the lower arm is the reference arm. The two arms enclose a magnetic flux <f>. The interferometer is 
connected to two transport leads, left (L) and right (R). 



Quantum Dot 




Reference arm 

FIG. 1: (Color online) Scheme of the Aharonov-Bohm intereferometer. 



We consider multiple electron transport channels, denoted with m. For every pair of incoming and outgoing channels 
(m,m'), there is a tunneling process through the dot that has a tunneling amplitude t m >_ m . The phases (f> m ',m = 
Arg(t rn i m ) are called intrinsic dot phases. For the tunneling amplitude through the reference arm, t re f = \t re f\e'" l ' r 'f , 
the phase can be set equal to the magnetic flux through the ring <p re f = $ (and in the numerical calculations we will 
take \t re f \ = 1). 

We are interested in calculating the total transmittance through the open interferometer, which is the experimen- 
tally available quantity. Incoherent leads are considered, for which the incident electron comes through each of the 
channels m = 1,...,M with equal probability 1/M. After a straightforward calculation (see Appendix A), the total 
transmittance through the interferometer can be expressed: 

Tint = Tdot + \Uef\ 2 + Jj\t re f \ \t m ,rn\ COs(0 m , m - <f>), (2) 

m 

with T dot = 1/AfEmro' \tm, m >\ 2 - Now we define the quantity t phase : 



. , y ] tm.rin (3) 



tphase — 

M 

m 

and its phase 4> e xp = Arg(t p hase)- With these notations, the total transmittance reads: 

Tint = T dot + \Uef\ 2 + ^\tref\\tphase\cOs{4> exp - $). (4) 

The above formula resembles Eq.l, but it is no longer the simple interference of just two waves. In the experiments 
0, El, HH , the transmittance through the interferometer exhibits AB oscillations when $ is varied and from the shift 
of the oscillations (when a gate is applied on the dot, for instance) one extracts the evolution of the measured phase. 
From Eq.4, it is clear that the measured phase is, in the multichannel case, <fi exp which can be expressed: 

^ ^-^m I COS 4>m,m /g\ 



VfEm \tm,m \ COS 4> m , m ) 2 + (J2,n l*m,m| Sin ( 



'm.m) 



The above formula is the main formal result of the paper. When the transport is single-channel, it is obvious that 
the measured phase 4> exp is equal to the intrinsic dot phase; for the multichannel transport, the above formula gives 
the value of the measured phase <j> exp in terms of the dot intrinsic phases § m ,m and amplitude modules |t m , TO |. 
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Eq. 5 already allows us to make a comment regarding the evolution of (j) eX p, if we assume that t m , m are of resonance 
type, with the module presenting a maxima at the resonant energy and the phase increasing, for instance, from to 
7r on the resonance width. Then between consecutive (in-phase) resonances, one has to add two complex numbers 
that are out of phase. This leads to a smaller amplitude of t p h ase , but also to a reduced variation of its phase, in the 
regions where the resonances overlap. We will show this explicitly in the example presented in the following. 



III. QUANTUM DOT WITH A MAGNETIC IMPURITY 



In this section and the next one, we will describe in detail an example of multichannel scattering, with the aim to 
obtain in the end the formula of the measurable phase (with the standard experimental method 0, El, EH), f° r the 
example considered. 

We consider here a quantum dot containing a S = 1/2 magnetic impurity , which interacts with the itinerant 
electrons via an exchange interaction. First, one must define the scattering channels and calculate the tunneling 
amplitudes. According to the standard scattering theory, the impurity dot represents a target that possesses internal 
degrees of freedom (the impurity may flip its spin) and, in such a case, the scattering channels (|m)) are the states of 
the electron+impurity system: | f-ff) , | |J|) , | fJJ.) , and | ) , where the simple arrow stands for the spin of the itinerant 
electron and the double arrow stands for the spin of the impurity (see also [25[ , where the same basis is used to study 
the spin transmittance through a closed interferometer). 

The Hamiltonian of the dot corresponds to the spin interaction between the electron and impurity: 



H D = -JiS=-J{s z S z + \s+S- + ]-s-S + ) 



(0) 



The terms in the Hamiltonian ca be written in detail: 

1 



s z s z = -\n){n\ + ^\u)in\ 



1 



7Ufr>U1H-7ltWtH 



s+s. = \n)iu 
ss + = \w(n 



(7) 



The quantum dot site index {no in Appendix B, where the general case is addressed) can be omitted, since the dot 
is single-site, and there is no risk for confusions (we write for instance | fit) f° r the dot site, instead of \hd, tH 1 ))- 



A 




A 



FIG. 2: (Color online) Scheme of the quantum dot containing a magnetic impurity, which interacts with the itinerant electrons 
through an exchange interaction J. The dot is connected to discrete leads. 



To calculate the tunneling amplitudes, we connect the dot to discrete leads, see Fig. 2 (for a comparison between 
discrete and continuous models in mesoscopic physics, see [1^|). The hopping constant between the dot and the leads 
is r, and the hopping between the leads neighboring sites i s 77 (in the numerical calculations we will take 77 = 1). 

Applying the general recipe in Appendix B (see also [13, [3] for equivalent demonstrations) for the particular case 
of a single site dot and four transport channels, one writes the effective Hamiltonian (which later shall be used to 
express the tunneling amplitudes): 





-J/4H 


-x 





" 


H eff = 





-J/4H 


-x 











J/4 + x 


-J/2 










-J/2 


J/4 + x. 



(8) 



where x = 2r 2 e %k and the columns and rows correspond to the following counting of channel vectors: | t1t)i| 14} > I t-H-) 
and I U). 
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Finally, as shown in Appendix B, the tunneling amplitudes between the channels m and m! can be expressed as: 



2i — sin(fc)(m'| — 

Tl E 



ttD 

U eff 



(9) 



There are six possible processes, corresponding to the non-zero elements of 1/(E — H®f), and respecting the total 
spin conservation. To simplify the index notation, the non-vanishing tunneling amplitudes will be denoted by U with 
i = 1 • • • 6: 



1) 
2) 
3) 
4) 
5) 
6) 

From the spin rotation symmetry of the dot Hamiltonian, one has: 



in) - 


Hit) 




= h = 


l*i 




in) - 


► in) 


tn,n 


= t 2 = 


1*2 




in) - 


► in) 




= t 3 = 


1*3 




in) - 


► in) 




= U = 


\u 


e i<pi , 


in) - 


► in) 




= *6 = 


1*5 




in) - 


► in 




= h = 


1*6 


e^ 6 . 



(10) 



tl = t(i, 
h = ti, 
*3 = *5- 



111) 



During the processes 3 and 5, the spins flip (see [29| for a discussion on the importance of spin- flip processes for 
open and closed interferometers), while the others are non-flip processes. 

In Fig. 3 we plot the amplitude modulus and phases for the three tunneling processes t\, t 2 and t 3 versus the 
energy of the incident electron. One notices that t\ is a single resonance process, while t 2 and t 3 have two resonances 
each. The single-resonance process, t\, has the classical Breit-Wigner structure, and the phase increases by 7r on the 
resonance width. 
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FIG. 3: The modulus and the phase (in 7r units) of the tunneling amplitude for the three processes a) ti, b) t-i and c) t$, 
plotted versus the energy of the incident electron. The parameters are: J — —0.5, r = 0.15. 

The processes with two resonances have different structures and they are of two types: with the two resonances 
in-phase (t 2 ) or out-of-phase (t 3 ). For the in-phase case, plotted in Fig. 3b (t 2 ), the phase increases from — tt/2 to 
approx 0.37T for the first resonance, then decreases between the resonances to the value — 0.37T and finally an increase 
to 7r/2 is noticed on the second resonance. In this case, the tunneling amplitude between resonances is lowered, 
approaching zero. For the out-of phase resonances, Fig. 3c (£3), the phase evolution is different, in the sense that the 
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phase evolves from tt/2 to 3ir/2 for the first resonance, and from 3tt/2 to 5ir/2 for the second resonance. Between 
resonances, the phase is roughly constant and the tunneling amplitude is increased. 

The above discussed behavior of the tunneling amplitudes can be understood if we perform a spectral decomposition 
of the dot Hamiltonian. It is straightforward to show that the dot eigenstates are: one singlet state with the energy 
Es = 3 J/4 and three triplet states with Et = — J/4. We define ts as the tunneling amplitude for the singlet state 
and tr the tunneling amplitude for the triplet states. One has ts = ir 2 sin(fc)((t-lj. | - \)(E — iT^L) -1 (| f-IJ-) - I it 
)) = 2ir 2 sm(k)/(E - 3J/4 - 2T 2 e~ lk ). Similarly, one obtains t T = 2ir 2 sm(k)/(E + J/4 - 2T 2 e~ ik ). The tunneling 
processes discussed above, £1,2,3, can be expressed in the terms of singlet and triplet (this results straightforwardly 
from the definition Eq.9; see also (30j): 

t\ = tT, 

t 2 =^(*r + *s) in phase, (12) 
t 3 = i(f T - ts) out of phase. 

The sign between the triplet and singlet amplitudes (tr and ts) in Eqs.12 explain the in-phase and out-of phase 
behavior seen in Fig. 3b and c. 

Now we can calculate the total transmittance of the dot, T^ot, by summing for all incoming channels (the probability 
of finding the incoming electron and the impurity in a given spin configuration is 1/4): 

T do t = \jZ\U\ 2 ^\jZ^\ 2 - ( 13 ) 

i=l " i=l 

Furthermore, by using Eqs.12, the dot transmittance becomes Tdot — l/4|ts| 2 + 3/4|£t| 2 , showing that the singlet 
resonance is three times lower in amplitude than the triplet resonance (see Fig. 4a, where for negative J the singlet 
resonance is lower in energy), which is an expected result for the two-spin scattering. We mention that this 1/3 ratio 
in the resonances amplitude was also found by Rejec et. al.[3l[ for the case of two electrons scattering and the result 
was used to give a possible explanation of the 0.25 and 0.75 anomalies in quantum wires. 



IV. THE MEASURED PHASE OF THE MAGNETIC IMPURITY DOT 



Let us now place the impurity dot in one of the arms of an open Aharonov-Bohm interferometer. The total 
interferometer transmittance is the one given by Eq.4, where t p hase = \{ti + ^2+^4 + t$) = \(t\ + £2)- Using the 
spectral decomposition (Eqs.12), one can write equivalently: 

tphase = ^(3- t T +t S ). (14) 

In Fig. 4a we plot the quantum dot transmittance Tdot, given by Eq.13, and the measured phase (f> exp = Arg(t p hase) 
versus the energy of the incident electron. One can notice that the phase has a reduced variation (by less than tt) 
both on the dot resonances and between them. 

The Aharonov-Bohm oscillations of the interferometer transmittance Ti nt are plotted in Fig. 4b versus the magnetic 
flux, in order to reproduce the experimental method @, [ll[ • Then, from the horizontal shift of the Aharonov-Bohm 
oscillations, one extracts the quantum dot "measured" phase evolution, plotted in Fig4a. The dotted vertical line 
in Fig. 4b is guide for the eye to see the phase evolution. For the lowest curve (corresponding to E = —0.6) one 
reads a phase 4> exp of about —tt/2 at zero flux (this means that we have neither a maximum of the oscillation, nor a 
minimum, but a value in the middle). The phase then increases, as the energy of the incoming electron increases, till 
the oscillations present a maximum at $ = zero (see the fifth curve from bottom); this means <j> exp = (from Eq.4 
one sees that, at zero magnetic flux, $ = 0, the oscillation presents a maximum if, and only if, <f> exp = 0). The phase 
continues to increase a bit more, and then the AB oscillations begin to shift in the opposite direction, meaning that 
the phase starts to decrease, and so on, resulting the phase evolution in Fig. 4a. For the considered parameters, the 
measured phase increases with ~ 0.657T on the singlet resonance, then decreases with ~ 0.557T between resonances, 
to finally increase a gain with ~ 0.97T on the triplet peak. Even if the shape of the phase evolution does not coincide 
with the features in 

evolves with less than tt is recovered by our multi-channel model. 



Ill ] - and a more realistic many-body model may be required - the fact that the measured phase 
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FIG. 4: a) The transmittance through the quantum dot with a S = 1/2 magnetic impurity (Tdot), and the measured phase 
(0ercp, given in in 7r units), plotted versus the energy of the incident electron; the parameters are: J = —0.5, r = 0.15. b) Total 
transmittance of the interferometer (Tint) versus the magnetic flux (in flux quanta) for the energies E G [—0.6, 0.4] marked with 
dots in (a) (E — —0.6 corresponds to the lowest curve). The curves are shifted on the vertical for the clarity of presentation. 
From the horizontal shift of the Aharonov-Bohm oscillations, one can extract the value of the measured phase for different 
energies. 



V. CONCLUSIONS 



An experimental method to extract the quantum dot phase , by monitoring the Aharonov-Bohm oscillations of an 
open interferometer, was developed in previous works (e.g. [3. [TH fl^|). 

In this paper wc show that, when the transport is multichannel, the phase extracted by using the mentioned method 
is not an intrinsic dot phase. Instead, it is the phase of the sum of the diagonal dot tunneling amplitudes. 

We payed a particular attention to the situation when an electron is scattered by a potential that mixes the 
transport channels. To this purpose, we have analyzed the particular case of a quantum dot containing a S = 1/2 
magnetic impurity. For this simple model, we have found variations of the measured phase with less than w, both on 
the resonances and between them. This reproduces some features of a recent experiment fill ] and is, in our case, a 
consequence of the multi-channel transport. 
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Appendix A: Total transmittance through a multichannel open interferometer with embedded quantum dot 

The measurable quantity in interference experiments is the total transmittance through the interferometer. To 
express this, we need the scattering solutions of the Schrodingcr equation. First, we chose the wave function in the 
left lead to describe the physical situation when an incident electron is coming from the left lead in a selected channel 
m and is reflected back with the amplitude r m i m into the channel ml: 

\^™) = J2e ikn \n,m}+ £ e- ikn r m , , m \n, m') . (Al) 

The wave function in the right lead describes the transmitted electron and for the open interferometer it is considered 
to be the sum of the two wave functions that correspond to the two paths: reference and QD arms. The tunneling 
amplitude through the reference arm does not depend on the transport channel m, t re f = \t re f\e l ^ r£f , and the 
corresponding wave function is: 

l*£/)=W£ e " ifc >> m >' (A2) 

On the contrary, the tunneling amplitudes for the electron that passes the quantum dot arm depend on the incoming 
and outgoing transport channels m and ml , and the wave function reads: 

1*0^)= E t m i, m e- ikn \n,m'), (A3) 

where t m i )OT is the transmission amplitude of the quantum dot for an electron coming in channel m and transmitted 
in channel m'. 

The interference wave function in the right lead of the AB interferometer will become: 

\* R n = i**7> + i^o 

= £(W +t m , m )e- lkn \n,m) + Yj t m ,, m e- lkn \n,m'). (A4) 

n£R n£R,m' ^m 

For the discussed situation, when the incoming channel in the left lead is m, the transmittance of the interferometer 
is noted by T m and is defined as: 

T m = \V R ' m \ 2 = \t ref + t mim \ 2 + K>, m \ 2 - (A5) 

Now we consider an incoherent lead for which the incident electron may come through each of the channels m with 
equal probability 1/M. The total transmittance becomes: 



Tint — j. i ^ ^ T m ^ j ^ "j [|^m,m'| ^m,m'\^ref\ 
rn m,m' 

)■ (A6) 



or, equivalently, 



2 

T int = T d ot + \Uef\ 2 + T^IWl £ Km, m| COs(0 m , m - $), (A7) 



where T dot = l/ME ro , m , \t m , m -\ 2 . 

The last term in the interferometer transmittance formula accounts for the quantum interference and depends only 
on the diagonal dot tunneling processes t m m that do not change the transport channel. The off-diagonal tunneling 
amplitudes t m ^ m i with m ^ ml add up in module in Eq.A7 and their phases do not influence the total transmittance 
(nor can they be measured). 
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Appendix B: Multi-channel transport through a quantum dot 

This Appendix aims to give an example (for the completeness of the paper) of how the tunneling amplitudes 
through a quantum dot can be calculated in terms of the dot Hamiltonian. See for instance [13, HH for a similar 
demonstration using the Green functions approach. The method presented here calculates the scattering solution for 
the system wave function. Although it is a rather general method (for the case when one electron is scattered by a 
region with interaction and internal degrees of freedom), it is practical only for few-body problems. For many-body 
problems, the task of calculating the system wave function is not realistic, and one needs different approaches. 

In this section, we use the same notation as in Appendix A and consider that the incident electron comes from the 
left lead through the channel m. The wave functions for the electron in the left and the right leads were already given 
in Eqs.Al,A3 and the wave function for the electron in the dot is: 

\V D - m )= E Cm>, m (n D )\n D ,m'). (Bl) 

n D £D dot ,m' 

ri£) denotes a site in the quantum dot. The coefficients r, t (from Eqs. Al,A3) and c (from Eq. Bl) are parametrically 
dependent of the incoming channel to, and will be determined form the Schrodingcr equation: 

Hd Hdl Hdr 
H LD H L 
Hrd H R 

The leads Hamiltonian correspond to the discrete chain with the hopping between nearest neighbors taken equal 
with the energy unit: 

H L{R)= Y, (\n,m'){n + l,m'\+h.c). (B3) 

neL(R),m' 

The coupling Hamiltonians are assumed not to change the channel, and account for the hopping of the electron 
between the dot site labeled nr> — ®l(r) an d the site of the left (right) lead that is nearest to the dot and is labeled 
1 (see Fig. 2) (like in Eqs.Al-3 the lead site index does not need to carry the L or R- subindcx, there is no danger of 
confusion) : 

H D ,L(R) = T D , L (R) Y \^L(R),m'){l,m'\, 
m' 

H L (r } ,d = T^ L(fi) E|l,?n')(0 L(i?) ,m'|. (B4) 

m' 

In the following, the coupling coefficients will be taken real and equal each other tdl = Tzji?. = t. 
From the solution in the lead, it is straightforward that the eigenenergy is E = 2cos(/c) (this is the dispersion 
solution for a discrete lead). 

In order to calculate the system wave function, we need to first write explicitly the set of three equation from Eq.B2. 
The steps are outlined below: 

1. One of the equations reads: Hld^ + HjJ^ L = E^ L . We use the explicit expression of wave function ^> L 
(Eq.Al) and write this equation for the site 1 of the left lead: 

{e 2lk \l,m)+Yl r ™\me~ 2lk | 1 , m'» + r ^ c m ,, m (0 £ ) 1 1, m) = 

m' in' 

E(e lk \l, m) + Y r m ', m e- ik \l, m')). (B5) 

Using Eq.A3 for the right lead we have: 

Ew, m e- 2lfe |lm')+TEc m -, m (0 fl )|l,TO') =£^t m( , m e-*|l,m'). (B6) 

m' m' m' 

One multiplies with (l,m'| on the left of Eqs. B5 and B6, and we have a relation between the dot wave function 
coefficients c m >. m and the lead wave function coefficients, r m > m and t m > m : 





\* D ) 








\* L ) 


= E 










1**) 



(B2) 
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I'm 1 ,m ^"^m / ,m(0L) ^m,m' ■ (■^'0 

2. The equation for the dot function reads: Ho^ D + tld^i + tbd^i = E^ D . We introduce the explicit form of 
* L and from Eq.Al,A3 and we have: 

Hd\*°) +T(e lk \0 L ,m) +Y / r ^, m e- lk \0 L ,m')) + rj] i^e"^, m'> = E\*»). (B8) 

m' m' 

Using the relation between the coefficients of the lead wave functions and dot wave function from Eqs.B7 one obtains 
(after some straightforward algebra) an equation for the coefficients c m < jm , which can be written in the compact form: 

(E-H% f )\*°) = 2iT(sink)\Q L ,m), (B9) 

with 

H° ff =H D + T 2 e- tk \0 a ,m')(0 a ,m'\ (BfO) 

m',a— L,i? 

3. Finally, the quantity of interest, the channel-dependent transmission coefficient, is: 

t m ,. m = 2iT 2 sm(k)(0 R ,m'\ 1 |0 L ,m). (Bll) 

h - H eff 
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